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Abstract. This paper treats the dynamics and scattering of a 
model of coupled oscillating systems, a finite dimensional one and 
a wave field on the half line. The coupling is realized producing 
the family of selfadjoint extensions of the suitably restricted self- 
adjoint operator describing the uncoupled dynamics. The spectral 
theory of the family is studied and the associated quadratic forms 
constructed. The dynamics turns out to be Hamiltonian and the 
Hamiltonian is described, including the case in which the finite di- 
mensional systems comprises nonlinear oscillators; in this case the 
dynamics is shown to exist as well. In the linear case the system is 
equivalent, on a dense subspace, to a wave equation on the half line 
with higher order boundary conditions, described by a differential 
polynomial pipx) explicitely related to the model parameters. In 
terms of such structure the Lax-Phillips scattering of the system 
is studied. In particular we determine the scattering operator, 
which turns out to be unitarily equivalent to the multiplication 
operator given by the rational function —p{iK)*/p{iK), the incom- 
ing and outgoing translation representations and the Lax-Phillips 
semigroup, which describes the evolution of the states which are 
neither incoming in the past nor outgoing in the future. 



1. Introduction 

In this paper we investigate the spectral theory, dynamics and Lax- 
PhiUips scattering for an abstract system which models the interaction 
between a finite dimensional linear subsystem and an infinite dimen- 
sional wave field on a halfiine. We will call such systems generalized 
Lamb models in that they extend the standard Lamb model (see [lOj ) 
that will be introduced shortly. Although our main concern is with 
linear oscillators, we will describe some properties of the models in the 
anharmonic case also. 

To introduce our models, let us consider a n-dimensional lagrangian 
system linearized around a certain equilibrium point. Its equations of 
motion are given by 

Gy = Hy 
1 
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where y G M" is the displacement from the given equihbrium point (for 
reference, y — 0), and the matrices G and H represent the quadratic 
approximation of kinetic and potential energy around the equilibrinm 
point. G is positive definite and both matrices are symmetric with 
respect to the standard inner product in M". For technical and theo- 
retical reasons it is more useful to endow M" with the inner product 
given by G. With respect to this inner product the matrix L — G^^H 
is symmetric and the lagrangian equation take the form 

y = Ly 

with L symmetric with respect to the G inner product. The case of 
a chain of harmonic oscillators is well known and yields to a Jacobi 
matrix for the operator L. 

Analogously, let us consider the wave equation on the halfline. To 
be definite let us consider Neumann boundary condition at the origin. 
Denoting with Ajv the 1-dimensional laplacian with homogeneous Neu- 
mann boundary conditions at the origin, the wave field (we have posed 
equal to one the wave velocity) evolves according to the wave equation 

(f) = Ajv0 . 

So we have two decoupled second order equations for two different 
oscillating systems, the finite dimensional one with generator L and 
the infinite dimensional one with generator Ajy. 

Thus on the direct sum L^(M_|_)©M" we have the self-adjoint operator 

Ao = Am®L 

and the corresponding abstract wave equation 

= . 

In a heuristic way a coupling between the two oscillating systems could 
be given by posing a constraint between boundary values of the wave 
field at the origin and the displacement of the finite dimensional sys- 
tem. The prototype of this coupling is fournished by the well known 
Lamb model where a semi-infinite string is coupled to a single particle 
oscillating in the transverse direction (see section 4.1 for the general 
case of a chain of oscillators); the particle, with mass M, is subjected 
to the tension T of the string at the origin and to a restoring harmonic 
force with spring constant K, so that the formal equations are given 
by the system 

^{t,x) ^(l)"{t,x) x>0, 
My{t)^-Ky{t)+T<P'{t,0+). 
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plus the constraint 

0(t,O+) = y(t). 

This model was proposed by Horace Lamb in 1900 as an example of 
dissipation in (subsystems of) conservative systems. In fact it is possi- 
ble to decouple field and particle dynamics, and the particle component 
satisfies a reduced equation which turns out to be, for t > 0, 

My{t) + 2Ty{t) + Ky{t) = T{<P'^{t, 0+) + 0o(t, 0+)) 

The forcing term on the right hand side depends on the evaluation at 
the origing of the free evolution for the wavefield of the initial data 
00, 00- Thus for initial data of compact support the forcing term is a 
pure transient definitively vanishing, and the reduced dynamics for the 
particle coincides for large times with that of a damped harmonic oscil- 
lator, so that the effect of interaction between particle and field reduces 
to damping only. This means exponentially fast return to equilibrium 
of the finite dimensional subsystem and correspondingly a neat transfer 
of energy to the field. This relaxation property towards the equilibrion 
position of the finite dimensional component is always true when the 
corresponding self-adjoint operator has empty point spectrum (see Re- 
mark EU]) as it is the case in the Lamb model. The result has various 
generalizations to anharmonic oscillators (see [8]). 
Some tridimensional models reduce themselves to generalized Lamb 
models due to simmetry. The case of an elastic spherical shell coupled 
to the acoustic field when radial oscillations only are allowed is treated 
in section 4.3, and it yelds to a nontrivial generalized Lamb model. 
Another issue of interest of these coupled systems are given by the fact 
that some linear models of classical and quantum field theory reduce 
themselves in the ultraviolet limit, and after due renormalizations, to 
generalized Lamb models. For example, the Schwabl-Thirring (see p6] ) 
model when restricted to its monopole sector (the only one where it 
is not trivial) and after a spring constant renormalization turns out 
to be equivalent to a Lamb model (see p^). A similar phenomenon 
occurs for the Pauli-Fierz model describing the interaction of a charged 
oscillator with the electromagnetic field in dipole approximation and 
after mass renormalization (see [3]) . In this case, however, reduction 
of the dynamics on its non trivial part, yields a boundary condition 
different from that of the Lamb model (see Example 14.21) . 

This discussion of motivating examples, and relevant studies existing 
in the literature, shows however that the usual formulation is partly 
formal in that it is not clear what it should be the functional setting 
of the Lamb system in the first place, and secondarily its Hamiltonian 
formulation, if any exists. 
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A guide to set rigorously these questions in this and in more general 
situations is suggested by an analysis of the coupling between field and 
particle. The idea is to restrict the uncoupled vector operator Aq to 
the linear variety defined by the constraint existing between field and 
subsystem; the uncoupled operator on this linear variety is no more 
self-adjoint but it is symmetric with defect indices (1,1). AH possi- 
ble selfadjoint extensions different from Aq itself correspond to a well 
defined coupling or interaction between two subsystems. The case of 
the Lamb model, for example, corresponds to the closed linear variety 
0(0+) = y. The most general boundary conditions still producing a self- 
adjoint operator, as we will see, is of the kind 9(f)'{0+) + 0(0+) = w ■ y, 
where ^ G M and w is a. given vector in M". The case of a chain of 
harmonic oscillators one of which coincides with the boundary point 
of the string, which is the more obvious generalization of the Lamb 
model, corresponds to a vector w with just a single entry nonvanishing 
and to 6* = 0. A generic w corresponds to nonlocal coupling between 
string and more than one oscillator, i.e. the interaction is not "nearest 
neighbour" . Our first concern, is to give a rigorous account of this con- 
struction and to explicitely describe the interacting system so obtained 
(see Theorem 12.11) as well as its spectral properties (see Theorem 12. 2p . 
The interacting operator so costructed is a singular perturbation of the 
selfadjoint operator Aq, related to the class of one dimensional point 
interactions, or better point interaction with inner structure previously 
studied in different context and with a different formalism by many au- 
thors (see e.g. [I3], [9], [1] and references therein). In passing, we note 
that the coupled operator we study corresponds to a boundary value 
problem for the wavefield only, but with an eigenparameter dependent 
boundary condition (see Remark 12.51) ; this sort of parameter depen- 
dent boundary value problems are well known in the literature, both 
physical and mathematical. However we do not follow this road to the 
study of spectral and scattering properties of the coupled operator. 

As a byproduct of the construction we obtain in Section 3 the Hamil- 
tonian structure of the system (see Theorem 13.21) . which we generalize 
to the case of anharmonic oscillators, giving conditions for the existence 
of global fiow (see Remark 13.31) . As far as we know, a completely rig- 
orous description of the Hamiltonian structure of such type of systems 
has been lacking up to now, whereas interesting, but formal treatments, 
are scattered in the hterature (see for example [12], [7]). 

In the case the symmetric operator L has no degenerate eigenvalues 
we show that the dynamics of the system is equivalent, for a dense 
set of smooth initial data, to a reduced dynamics of a wave equation 
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on the half line which incorporates the interaction with the finite di- 
mensional systems through a higher order boundary condition of the 
kind p{dx)4>{t, 0_|_) = 0, where the polynomial p is explicitely related 
to the paramenters entering into the definition of the model (see The- 
orem 15.11) . This is a technical result, useful for the analysis of the 
Lax-Phillips scattering for the system, which is the main topic of the 
remaining part of the paper. 

In Section 6, in the case of empty point specrum, we determine the 
incoming {R~) and outgoing (R^) translation representations which 
make the dynamics unitarily equivalent to the translation on L^(R) 
defined by T^f{x) := f{x — t). This provides the scattering operator 
S* for the system by the relation S* = R^{R~)~^. Moreover S* turns 
out to be unitarily equivalent to the multiplication operator given by 
the rational function —p{iK)*/p{iK). In Section 7 the Lax-Phillips semi- 
group Z*, t > 0, which describes the evolution of the states which are 
neither incoming in the past nor outgoing in the future is completely 
characterized. It acts on a finite dimensional vector space, whose di- 
mension coincides with the degree of the polynomial p, by = e~*^, 
where the spectrum of the generator B is made of the resonances of 
the system. Such resonances correspond to the roots of the polynomial 
P- 

2. Singular perturbations of the free dynamics 

Let us begin with some definitions. We denote by L^(M_|_) the Hilbert 
space of square-integrable functions on the half-line (0, +oo) and by 
i/^(M+) and H'^(R+) the Sobolev spaces 

H\R+) := {(f) G L^{R+) : 0' G L\R+)} , 
H\R+) := {0 G L\R+) : 0', 0" G L\R+)} . 

Here the prime 0' denotes a spatial derivative. With a dot, 0, we will 
denote a time derivative. We then define if|r(R+) as the subspace of 
H'^{R^) of functions which satisfy homogeneous Neumann boundary 
conditions at zero, i.e. 

HU^^) := {0 G H'{R+) : 0'(O+) = O} . 

We denote by (■, ■)2 and by || ■ II2 the usual scalar product and the 
corresponding norm on //^(M^). 

Given the n-dimensional Hilbert space f) with scalar product (■, ■) and 
corresponding norm || • ||, and given the symmetric operator L : f) — > {), 
we consider the self-adjoint operator 

Ao : HU^+) © i) C L\R+) © [) ^ L\R+) © i) , Ao(0, y) := (0", Ly) . 
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Regarding the spectrum of Aq one has 

o'essiAo) = aaciAo) = (-00, 0] , app{Ao) = a{L) . 

In order to couple the two dynanical sistems described by the equations 
(p = 0" and y = Ly we define the continuous and surjective hnear 
operator 

r : H\R^) © ^ C , r(0, y) := 0(0+) - {w,y) , w e i) , 

and then we consider the closed symmetric operator Aq obtained by 
restricting Aq to the kernel of r. Aq has deficiency indices (1, 1) and 
we are interested in its self-adjoint extensions different from Aq itself, 
which we parametrize by the real extension parameter 6. Thus to each 
quadruple (l),L,iy, 6*) corresponds a different generalized Lamb model. 
The next theorem completely characterizes such models. 

Theorem 2.1. For any 6 the linear operator 

A : D{A) C L\R+) © f) ^ L\R+) © [) 

defined by 

D{A) := {(0,y) G H\R^) © [) : ^0'(O+) + 0(0+) = {w,y)] , 

A(0,y) := {<P",Ly + w^'iO+)) 
is a self-adjoint extension of Aq and its resolvent is given by 
{-A + z)-' = {-Ao + z)-' + {9 + T{z))-'G, © G,* , 



where 



T{z) ■=-(^±^ + {w, (-L + z)-^w)^ , ±Rev^ > 



and 

G,= {±^-^, -{-L + z)-^w] , ±Rev^>0. 



Proof. We will make use of the mathematical procedure developed in 
[T4j (see also j4j, Theorem 2.2, for a similar proof in the case of a 
one-dimensional model in acoustics). 

For any z G p{Aq), let us consider the two linear continuous operators 

Giz) : L2(R+) © ^ C , Giz) := t{-Ao + z^^ , 
Giz) -.C^ L\R+)®^, Giz) := G{z*y . 



Since 



+ z] ■.L'^(R+)^H^ 
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has kernel 

gN{z;xi,X2) = ± , ±Rev/i>0, 

the operators G{z) and G{z) are represented by the vectors G^* and 
Gz respectivelty, where 

G-,= [± — -{-L + z)-^w\ , ±Rev^>0. 



Note that 

(2.1) Ran(G(^)) n D(Ao) = {0} . 
Now we define, for any z G p{Aq), the map 

r(z):p(Ao)-.C, T{z):=-tG{z), 

i.e. 

V{z) := - + (w, (-L + z)-^w)^ , ±Re^/^ > , . 

At first let us note that the the function F satisfies the relation 

(2.2) r{z)-r{w) = {z-w)G{w)G{z). 
Indeed 

T{z)-T{w) = t{G{w)-G{z)) 
and, by first resolvent identity and by the definition of G{z), 

{z - w) (-Ao + z)-^G{z) = G{w) - G{z) . 
Relation (12.21) implies that 

R{z) := i-Ao + z)-^ + {9 + T{z))-^Gz ® G,* 
satisfies the first resolvent equation 

(2.3) R{w) - R{z) = {z-w) R{w)R{z) . 

By the definitions of G{z) and G{z), and since T{z)* = T{z*), one 
obtains 

(2.4) R{z)* = R{z*) . 

Moreover, by (12. ip . the linear operator R{z) is injective. Thus 

A := -R{z)-^ + z 
is well defined on the domain 

D{A):= Range(i?(^)) . 
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By (12.31) such a definition is ^-independent. By (12. 4p A is symmetric 
and is self-adjoint since 

Range(-A ±i) = L'^{R+) © f) 

by construction. We have thus defined the self-adjoint operator 

D{A) := y,) + {9 + T{z))-\MO+) - {w, y.)) , 0, G , 

{-A + z){(f),y) := {-Ao + z){^,,y,). 

This implies 

0'(O+) = -{e + T{z))-\MO+) - {w,y.)) 

and 

0(0+)= 0,(0+) T^0'(O+). 

Therefore 

^0'(o+) = {6 + v{z)) 0'(o+) - r(^) 0'(o+) 

- 0,(0+) + {w, y,) +{^ + (-^ + \ 0'(O+) 



= - 0'(O+) + {w, {y, + (-L + zY^w) 0'(O+)) 
= -0(O+) + (w,l/). 
Posing 

A{(t),y) = (Ai(0,|/),A2(0,y)) 

one obtains 

A(0,i/)(x) = 0:(x)T#'(o+) 



2; 



(x)T0'(O+)^) =0"(x) 



and 

^2(0, y) = Ly, + #'(0+) (-L + 2;)-iu; 
=Ly + (-L(-L + z)-i + z{-L + 2)-^) 0'(O+) 
=L|/ + w0'(O+). 

Let us define the (eventually empty) set 

or^(L) := {A G a(L) : w G [)f } , 
where f)A denoted the spectral subspace relative to A. 



□ 
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Theorem 2.2. 

aess{A) = (TaciA) = (-00,0] , 

app{A) = a^{L) U {A G p{L) n (0, +oo) : 6 + r(A) = 0} . 

Proof. The properties regarding the essential and continuous spectrum 
are more or less standard and can be obtained proceding as in [[4], 
Theorem 2.3. Let us now come to the point spectrum. 

1. Let A e cr(L). Then (0, y\) is an eigenvector if yx solves the equations 

Lyx = Xyx, {w,yx)=0, 

thus A G (TwiL). 

2. Let A > 0. Then (f)x{x) := e"^^' solves (pl = X(f)x- Thus {(j)x,yx) is 
an eigenvector if A and yx solve the equations 

(2.5) Lyx- y/\w = Xyx, 

-eVX + l-{w,yx) = 0. 

If A G p{L) then 

yx = -^/X{-L + X)-'w 
and A must solve the equation 

-e\^ +l + VX{w, {-L + X)-^w) = . 

If otherwise A G <j{L) then (12.51) can be solved only ii w & i)-^ by 
yx = yl + yi, where ylei)x and G [)f is defined by 

yi := - i-Lx + Xr'w , Lx := (1 - Px)L{l - Px) : i)i ^ i)i ■ 

Thus A G o"«,(-L) and moreover it has to solve the equation 

-e^/x + l + VX{w, {-Lx + X)-^w) = . 

□ 

Remark 2.3. When aw{L) is empty, i.e. in the generic situation, the 
point spectrum of the interacting operator A is quite different from 
the point spectrum of the decoupled one, Aq. In particular, the free 
eigenvalues of the finite dimensional subsystem disappear, and in their 
place could possibly appear the real solutions of the equation r(A)+6' = 
0. In fact, as we shall see in Section 6, the disappeared eigenvalues, 
which for the uninteracting operator Aq are immersed in the continuum 
spectrum, become resonances of the interacting operator. 
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Remark 2.4. In the case (Jw{L) 7^ we can suppose, without loss of 
generahty, \) = \)i ® L = Li ® L2 and w = Wi ® Then for the 
self-adjoint extensions given in Theorem 12.11 we have A = Ai (B L2, 
where 

Ai : D{Ai) C L\R+) © f^i L\R+) © f)i 

is defined by 

D{A^) :={(0,yi)G//2(M+)©f)i : e0'(O+) + 0(0+) = yi)^ J , 

(0,1/1) := (0",Liyi + w7i0'(O+)), 

i.e. the dynamics on 1)2 is trivial, in the sense that it is decoupled 
from the field one. Thus the hypothesis crw{L) 7^ is equivalent to the 
existence of a subspace on which the particles and the field are uncou- 
pled. In other words the points of the pure point spectrum belonging to 
(jyj^L) correspond to " radiationless motions", in which the interaction 
between oscillators and field is decoupled. Similar exceptional solu- 
tions in which the subsystem oscillates at a normal frequency of the 
decoupled system are known, for example, also in the classical electro- 
dynamics of an extended charge where they are called Bohm-Weinstein 
modes. In that case the coupling between field and particle is defined 
by the charge density p{x) of the particle, and the condition to have 
radiationless modes of frequency uj is that the Fourier transform of the 
form factor satisfies f){uj) = 0. 

Remark 2.5. The operator A can be interpreted in a formal way as a 
differential operator with an eigenvalue dependent boundary condition. 
Let us consider the secular equation for the operator A and in particu- 
lar its finite dimensional component, and couple it with the boundary 
condition for elements of the domain of the operator. We get 

Ly + w(f)'{0+) =Xy 

From the first equation it follows y = — 0'(O+)(L — X)^^w and substi- 
tuting in the second equation one gets 

{9 + {w, {L - \)-'w)) 0'(O+) + 0(0+) = 

which is, formally, a Robin boundary condition for the field at the 
origin. The condition contains the eigenvalue A and it is is known in the 
physical and mathematical literature as an energy dependent boundary 
condition. From this point of view, the boundary value problem for the 
coupled operator can be reduced to a boundary value problem for the 
field only, but eigenvalue dependent. 
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By the way note that, as it should be, the above eigenvalue dependent 
boundary condition is equivalent to the eigenvalue equation in Theorem 
2.3, r(A) + ^ = 0, as it is immediately seen by the position = A. 

3. The Hamiltonian structure 

In this section we are interested in describing the Hamiltonian struc- 
ture of the dynamical system related to the abstract wave equation 



The solution of the corresponding Cauchy problem is then given by the 
symplectic flow generated by the linear operator [. Y\. 



We begin by determinig the quadratic form corresponding to —A: 
Theorem 3.1. Let us denote by Q the quadratic form of —A. 



1. If 6 = then 

DiQ) = {{cj>,y)eH\R+)®i^ : cP{Q+) = {w,y)} , 

Q:D{Q)^R, Qi^,y) = mi-{Ly,y). 

2. Ife^O then D{Q) = H^{R^) ® I) and 

Q : H\R+)®^ ^ M, Q(0,|/) = y\\l-{Ly,y)-\ |0(O+) - {w,y)\' . 



Thus the proof is done if Q is bounded from below and closed. This 
follows from 





{w,y)) 




a>0, 0<6<1. 



□ 




and deflne 



Ho := D{Q) ® L2(R+) ® ^ . 



Then one has the following 



12 



MASSIMO BERTINI, DIEGO NOJA, AND ANDREA POSILICANO 



Theorem 3.2. The linear operator 
1 



A 



: D{A) ® D{Q) C Ho ^ Ho 



is the generator of a strongly continuous group of evolution 

f/* : Tio — ^ Tio 

which preserves the energy 

E{{<p,y),{^,y)) :=i (q(0, ^) + ||0||2 + 

Such an operator is the Hamiltonian linear vector field corresponding 
to the quadratic Hamiltonian E via the canonical symplectic form on 
L\R+) © i) © L\R+) © i) given by 

^((01, yi, 01, yi), (02, 2/2, 02, 2/2)) 

:=(0i, 02)2 - (02, 02)2 + (2/1, h) - (2/2, 2/2) • 

Proof. The operator A is self-adjoint and bounded from above. Thus 
the resuh concerning evolution generation follows from the theory of 
abstract wave equations (see e.g. [B], chapter 2, section 7). The re- 
sults about the Hamiltonian structure follows from the theory of linear 
Hamiltonian systems in infinite dimensions (see e.g. [5], chapter 2). □ 

Remark 3.3. The above results can be immediately extended to a 
nonlinear situation. Indeed, given the potential function V let us con- 
sider the Hamiltonian 

H : 7io — ^ M , 

where 

H{{<P,y), (0,2/)) = I {ml + ml + \\yr) + V{y) 
when 6 = and 

^((0,2/), (0,2/)) = I (\\<P'\\l + 11011^ + \\yr-l 10(0+) - {w,y)fyV{y) 
when 6^0. 

The non linear Hamiltonian vector field corresponding, via the canon- 
ical symplectic form on L^(]R+) © f) © L^(]R_|_) © f), to if is given by 

Xh : DiA) © DiQ) C Ho ^ Ho , 

Xh((0, y), (0, y)) := (0, y, 0", -VV{y) + w 0'(O+)) . 

Obviously Xh = Xe + B, where Xe is the linear Hamiltonian vector 
field corresponding to the quadratic Hamiltonian E and B is vector 
field B{{(f),y),{<i),y)) := {0,0, -{Ly + VV{y)),0). Thus if V is twice 
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differentiable then by Segal's existence theorem (see [TTj) Xh generates 
a local continuous non linear symplectic flow on TCo- Since Q is bounded 
from below, if 

V{y)>ci\\yf-C2, ci>0, C2 > , 
then such a flow is global. 

4. Examples 

Example 4.1. The dynamics of the Lamb model (see [lO]) we already 
described in the introduction, given by the equations 

ip{t,x)=<j)"{t,x) 
My{t) = -Ky{t)+T^'{t,0+) 
y(t)=</)(t,0+) 

is described by the self-adjoint extension A corresponding to 

M K T 

dim{) = l, (x, y) = — X*?/ , Ly = -—y, ^ = ^ 9 = 0. 

The similar model with n point masses 

0(t,x) = 0"(t,x) 

M,y,{t) = -K,{y,{t) - y^it)) + T0'(t, 0+) 
M2y{t)2 = -K2{y2{t) - ysit)) + K,{y,{t) - y2{t)) 

Mn-iyn-l{t) = -Kn^l{yn-l{t) - yn{t)) + (yn-2 (t) " l/n-l(t)) 

Mnjjnit) = -Knyn{t) + i^n-l (Z/n-l (i) " 2/n(t)) 
yi(t) = 0(t,O+) 

is described by the self-adjoint extension A corresponding to 

1 " 

dim f) = , (x, ?/) = - ^ Mj x*yj , 



/ Ml Ml ^ ^ ■ ■ ■ ^ 

-ft'i i^i+i<:2 ^ n n 

M2 M2 M2 ^ • • • ^ 

Q X2. x2+i^3 n 

^ Ms Ms Ms • ■ ■ ^ 



n J^ra-3 JVn-3-t-J^n-2 i^n-2 r\ 

^ ... ji/f (1^ ^ 





V 







M„_2 


M„_2 





A'n-2 


M„_i 









M„_2 

Kn-2+Kn-\ -ffn-1 



M„_i M„_i 

i^n-l Kn-l+Kn 

Mn M„ J 
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W= (^,0,...,0l , ^^ = 0. 



The corresponding Hamiltonian is given by 

H : {{(j),y) e H\R+)®C'' : (P{0+) = yi] ® L\R+) ® C" 

1 



Hii<f>,y),i<j,,y)): ^ 




1 " 1 \ 

k=l / 



where the matrix A is given by 

/ Ki -Ki ... 

-Ki K, + K2 -K2 ... 

A= : : 

... —Kn-2 Kn-2 + Kn-l — -R'n-l 

V ... -Kn^i Kn^i + Kj 

and ■ denotes the standard inner product in C". 

In the following examples we describe models, from classical elec- 
trodynamics and theoretical acoustic respectively, which are not inter- 
pretable as standard Lamb models in that 6' 7^ 0. 

Example 4.2. The renormalized Pauli-Fierz model. 

A three dimensional charged oscillator characterized by frequency uj, 
mass m and electric charge e interacting with the electromagnetic field 
in dipole approximation has a dynamics described, in the point limit 
and after mass renormalization, by a well defined self-adjoint operator 
which couples particle momentum and vector electromagnetic poten- 
tial. Dynamics and its main properties, classical and quantum, are 
constructed and studied in [3] . This is the point limit of the Pauli-Fierz 
model for the case of a quadratic potential energy. Due to the dipole 
approximation, the action of this operator is non trivial (i.e. different 
from the free uncoupled dynamics) only on the radial component of the 
field, and by standard decomposition using vector spherical harmonics 
it turns out that on this monopole subspace, the restricted dynamics 
for every couple {4>,p) constituted by a component of the vector po- 
tential on the non trivial subspace, and a corresponding component of 
the particle momentum, is given by the coupled system 

4>{t,r) = <f)"{t,r) 

p{t) = -—cu^^it,0+) 

0'(t,O+) + 00(t,O+) = ^p(t). 
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Writing the field which appears in the evolution equation for p in 
terms of its derivative 0' and p by means of the boundary condition 
one obtains 



(f){t,r 
Pit)-- 
2e2 



-4>"{t,r) 

-uj^p{t)+eLu'^(f)'{t,0+ 
2e 



3m 



0'(t,O+)+0(t,O. 



3m 



Pit)- 



Let us remark here that by Newton's law p = —mu'^q, where g is a 
component of the particle position. Thus the Cauchy initial datum for 
p is obtained from the initial position. 

In conclusion the dynamics of the renormalized Pauli-Fierz model in 
dipole approximation and with quadratic external potential is described 
by the self-adjoint operator A corresponding to 



dim f) = 1 , {x,y) 



2x*y 



Ly = -cj^ y , 



The corresponding Hamiltonian is given by 



w 



9 



2e2 
3m 



+ 



3m 



3m 



p 



As recalled in the introduction, a field-particle interaction which re- 
duces to the standard {9 = 0) Lamb model in the point limit and 
after spring constant renormalization is the Schwabl-Thirring model, 
in which a scalar field interacts with a scalar oscillator (for details, in 
a different framework, see [12]). 

Example 4.3. A spherical elastic shell in the acoustic field. 

Let us consider the exterior problem for a spherical shell of mass 
M, radius Rq and constant surface density p = undergoing radial 
motion only, and interacting with an irrotational acoustic field in the 
linear approximation. The shell is elastic, i.e. on every surface element 
acts a restoring force proportional to the radius, and of Young modulus 
K. If small radial oscillations around Rq are considered, introducing 
the variables -ip, related to the acoustic potential by (f){Ro + r) = 
r > 0, and the radius Rq + R{t) and taking into account the continuity 
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of velocity 0' of the acoustic field at the boundary one obtains the 
equations of motion (assuming propagation velocity equal to one) 

i'it,r) = i>"{t,r) 

MR{t) = -KR{t) +4:7iRoptp{t,0+) 

R R2— = ■ 

This system can be converted to a generalized Lamb system by intro- 
ducing the new variable (in fact a sort of total momentum) 

P := MR-4:nRopip{0+). 

Re-writing the above equations in terms of the new variable P, ex- 
pressing the field ip which appears in the evolution equation for P in 
terms of its derivative ip' and P by means of the boundary condition, 
one obtains 

V^(t,r) =<(t,r) 



^'{t, 0+) + 4^{t, 0+) = ,7 P{t) 



As regards the initial datum for P it can be recovered from the one for 
R as in the previous example. 

By defining u;^ := ^ the above system is described by the self-adjoint 
operator A corresponding to 

X* y up' 
diml) = l, (x,?/) = — — , Ly 



Anu^Rl ^ Ro 
w = — , — 







1 + Ro ' 1 + Ro 

The corresponding Hamiltonian is given by 

H : HUR^) © C © L^(R^) © C ^ 

/7((^,P),(^,P)): 
1 

An analysis of interaction of elastic surfaces with acoustic fields from 
the point of view of Lax-Phillips scattering theory is given in [2j. For 
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a study of one dimensional models in acoustics in the framework of 
self-adjoint extensions we refer to [1]. 



5. Wave equations with high-order boundary conditions 

From now on we will consider self-adjoint operators A which are self- 
adjoint extensions corresponding (according to Theorem 12 .11) to L's and 
w's such that 



(5.1) 



is a basis in I). 



Note that the examples given in Section 4 satisfy such hypothesis. 

With respect to the orthonormal base obtained from {L''w}q~^ by 
the Schmidt orthogonalization procedure the linear operator L is rep- 
resented by a Jacobi matrix. However we prefer to consider here the 
unitary isomorphism f) ^ induced by the orthonormal system {cj}" 
made of the eigenvectors of L. For any vector y & i) and for any linear 
operator M : f) — > [) we pose 



y = {yi, ■■■,yn) 



With these notations 



M 




w = {wi, 



where cr(L) = {Ai, 



L = 







\0 





A2 







XnJ 



A„}. We introduce the diagonal matrix 
.. 

W2 ■ ■ 



W = 











Wn) 



the Vandermonde matrix 



V 



( 1 

Ai 



1 

A2 



A 



n-l 



1 \ 



and then we define 



M := VW . 
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Since 

det M = J] (A,- - A,) J] Wi^Q, 

our hypothesis (15.11) is equivalent to 

Aj 7^ \j and cr^(-L) = 0. 

Thus under our hypothesis the spectrum of L is simple and, by Theorem 
12.21 A has no eigenvalue immersed in the continuos spectrum. 

Let us denote by iS(]R+) the space of rapidly decreasing smooth func- 
tions on [0, +oo). We define the dense subspace V C Ho hj 

V:={{(P,y,<P,y)eno : G 5(R+), G 5(M+), 

where 

n 

= ^Pkidx)(l) (0+) Cfc , 
k=l 

and Pk{dx) is the differential operator with constant coefficients asso- 
ciated to the polynomial recursively defined by 

pi{z) = ez + l, pk{z) = z'^pk-i{z) ~ {w,L''-'^w)z , k>2. 

The next theorem is the main technical point as regards the successive 
study of the Lax-Phillips scattering of generalized Lamb models. It says 
that V is invariant under the flow and that on this dense subspace 
a generalized Lamb model is equivalent to a wave equation with a high 
order boundary condition at zero. 

Theorem 5.1. Let the strongly continuous group of evolution pro- 
vided by Theorem \3.<X Then 

and 

= (0(t), M-'v{<P{t)), 0(t), M~'vi^m 
where (j){t, x) solves the equations 

dtt(pit, x) = dxx(pit, x), x>0, 
(5.2) p(a,)0(t,O+) = 

0(O,x) = 0o(a;), 0(O,x) = 0o(x) . 
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Here p{dx) denotes the constant coefficients differential operator of de- 
gree 2n + 1 (2n if 6 = 0) associated to the polynomial 

n 

p{z)=Pr.^,{z)-Y,K {V-'),^P,{z). 
Proof. Let be the solution of the Cauchy problem 

^,m,y{t)) = Am,y{t)) 

(^•3) (0(O),i/(O)) = (0o,l/o), 

(0(O),y(O)) = (0o,2/o) 

with (00, Vo, 00; Vo) £ ^ and let us suppose that (0(t), 0(i)) is in iS(M+) 
for all times. Then deriving with respect to time the boundary condi- 
tion 2fc-times, k = 0, . . . ,n, using 0(t, 0+) = 4>"(t, 0+), one obtains the 
n + 1 equations 

{w, L^y{t)) = Pfc+i(9,)0(t, 0+) , A; = 0, 1, . . . n . 

The first n of such equations can be can be rewritten as 

My{t) = vi^it)) 

so that 

y{t) = M-\{<P{t)) , y{t) = M-\{^{t)) 

for all times. Moreover, inserting the expression for y{t) into the n'th 
equation one obtains p{dx)(t)if) 0+) = 0, so satisfies f l5.2p . 
Conversely let 0(t) be the solution of (15.21) and put 

y{t) := M-^v{<P{t)) . 

Then the n equations 

{w, L^y{t)) = Pk+i{d.)(t){f 0+) , A; = 0, 1, . . . n - 1 . 

are satisfied. The first equation says that (0, y) and (0, y) are in 
D{A). Deriving each equation two times with respect to time and 
using 0(t, 0+) = 0"(t,O+) one obtains 

{w.L^yit)) =R-+i(5x.)0"(i,O+) =pfc+2(5.)0(t,O+) + (w, L'^w)0'(t, 0+) 
= (w, L^{Ly{t) + w(j)\t, 0+)) , A; = 0, 1, . . . n - 1 , 

which implies y{t) = Ly(t) +w0'(t,O+). So {4>,y) is the solution of 
(15.31) . This also show, by unicity, that if the initial conditions of (15.31) 
are in iS(]R4.) then they are in iS(M-|_) for all times. This justifies the 
assumption we made at the beginning of the proof. □ 
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The next lemma makes more explicit the polynomial p appearing in 
the previous theorem: 

Lemma 5.2. 

n / j \ 

p{z) = {Oz + 1) det(z' ("j-ki^^ ^^~^'^) -^"^""'^ ' 

j=i \k=i J 

where 

oo = 1 , ttj := (-1)^ ^ \i^...Xi., l<j<n. 

ii<---<ij 

Proof. Put 

n 
i=l 

and 

n 

bjk := - XI \^~^~Vir > l<A;<n+l,l<j<A;-l. 
1=1 

By the definitions oip{z), Pk{z), dj and bjk one has 

Ti+l / ^ 

p{z) = {ze + 1)Y: ~a,z^^'-'^ + W E Wi^^^^-^^ + E % E ^^^--'^'"^ 

n+1 n / n+1 \ 

= {ze + 1) E + ^ E E ^^k^k 

j=l j=l \k=j+l J 

^{ze + l)pa{z'')+Zpbiz^) , 

where 

n+1 n 

Po.^^) = E = E ' 

j=l j=0 

n n 



n n+1 n+1 

= - E E >^k''"~<^k? = - E ^^(^' ^'"'"'^) ■ 

fc=l i=j+l i=J+l 

Since 

n+1 n n 

E = - E E -^^^ (^"')^. ^^•'^ = K-K = ^. 

j=l j=l i=l 
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the eigenvalues Ai, . . . A„ are the roots of pa- Thus 

Yl i-^r~'^'K---Kn-,+., l<J<n, an+i = l- 

il<---<in-j+l 

□ 

The next result gives the relation between the roots of p and eigen- 
values and resonances of the self-adjoint extension A: 

Lemma 5.3. Suppose det L ^ and let us define the couple (0, y) by 

<^(a;) := e^"" , y := M'^v{(p) . 

Then 

p[z) = <^ I Ly + u;0'(O+) = z'^y ^ ze C\iR . 

[e^'(0+) + ^iO+) = {w,y), 

Hence 

p{z)=0, Re{z)<0, <^ z = -^/X, A G o-pp(A) n (0, +oo) . 

Proof. Since 4>" = z^4> implies 2; 7^ and p(0) = (— l)"detL 7^ 0, we 
can take z 0. By the definition of (0, y) we only need to show that 

p{z) = <^=^ Ly + wz = z^y'^ . 

Moreover, beside {w,y) = pi{z), one has 

{w, L'^y) = Pk+i{z) = z^pkiz) - {w, L'^^'^w) z, k = 1, . . . ,n - 1 . 

The above equalities, together with p{z) = 0, gives 

(W, L"7/) = Pn+l{z) = Z^Pn{z) - (w, L^'^w) Z 

Thus 

(w, L''{Ly + WZ- z^y)) = 0, A; = 0,...,n-1, 

i.e. Ly + wz — z^y. Reversing the above argument one has that 
Ly + wz — z'^y implies p{z) — 0. 

Suppose now piiv) = 0, G R, so that Ly -\- iuw — —v^y. Since 
(w, = iv9 + 1 we have 

{{L + u^)y,y)^u^e-w. 

Since L -|- i/^ is symmetric we have u — Q. But i/ 7^ by det L 7^ 0. □ 

Remcirk 5.4. By the previuos lemma we have that the polynomial p 
has no purely immaginary roots. The ones in the left half-plane are 
real and correspond to eigenvalues, the ones in the right half plane 
give rise to not normalizable solutions of the eigenvalue equations and 
correspond to resonant states. 
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6. Lax-Phillips scattering 
From now on we suppose that 
(6.1) app(A) = 0. 

Here we remark that the successive results hold, with the appropriate 
modifications, even without this hypothesis which is just a convenient 
choice in order to simplify the exposition. 

Since we already supposed crw{L) = 0, by Theorem 12.21 the above 
hypothesis means 

{A G p{L) n (0, +oo) : e + r(A) = 0} = , 

i.e. we are supposing that there are no strictly positive solutions x of 
the equation 

I " U„.|2 



X ' ^ —A,- + x^ 



9. 



This is true if and only if 

C (-00,0), e<0. 

Thus hypothesis (16.11) is satisfied by both Examples 14.11 and 14. 3[ 

By Lemma 15.31 (16. ip implies that the polynomial p has no nega- 
tive real root and that the complex ones (which appear in complex- 
conjugate pairs since p has real coefficients) are all contained in the 
right half-plane. Thus 

app{A) = Roots(p) C {zeC : Re{z) > 0} . 

By functional calculus, since A is injective and negative by (16.11) we 
have that 



cos 



-v^^sinty^ cost^/^. 
Moreover f/* extends to a strongly continuous unitary group 

where Ti is the Hilbert space given by the completion of Tio with respect 
to the scalalr product corresponding to energy norm 

||(0,y,0,y)|U:=E(0,l/,0,y)^/^ 

Remark 6.1. By Theorem 12.21 our hypothesis app{A) = says that A 
has purely absolutely continuous spectrum. Thus 

lim \\y{t) II = . 

t— >itoo 



LAX-PHILLIPS SCATTERING FOR GENERALIZED LAMB MODELS 23 

Indeed by functional calculus and Riemann-Lebesgue lemma, denoting 
by P{dX) the projection- valued measure corresponding to A, one has, 
for any {(l),y,^,y) e Ho, 

lim Xiyi{t) 

t—>±oo 



poo 

= ^hm^ / costy/X{V^{0,ei),P{dX)V^{(f),y))mR+)m 

J 

+ ^hm^ / siniVA(y=A(O,ei),P(dA)(0,y))i2(K+)ei, = 

J 



Analogously one has 



lim ||y(i)|| = 0. 



Define 

T {(/-,/+) e 5o(M) X 5o(M) : = 0} , 

where 

5o(M) :=|/e5(M) : / /(x)dx = 

^{/g5(M) : / = ^?', ^?g5(M)} . 

By considering the Fourier transform (denoted by " ) of the differential 
equation 

(6.2) p(a,)/_+p(-5,)/+ = 

one obtains 

Thus one has the following result, which permits to define what will be 
the scattering operator. 

Lemma 6.2. 

(/-,/+) ^ U{n) = -^^UK). 



Since p has real coefficients, wc have 

p(m) 



p{—iK) 



1, 



therefore 



S, : 5o(K) ^ 5o(M) , ( VH^) -^^^ /(/^) . 
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extends to an unitary map on L^(]R). 

Let (j)(t) be the solution of (15.21) with initial data (0, y, 0, y) G V. 
Then 

(j)(t, x) = a{x + t) + b(t — x) , 

where the couple (a, b) is determined on an half-line up to a constant 
c by 

r+oo 

(6.3) a{x) = --J {4>{y) + ct>'{y))dy + c, x>0, 

(6.4) h{-x) = -j (0(2/) - 0'(2/)) rfy - c , X > . 

The functions a{x) and b{—x) are then determined for the remaining 
values of X < by solving the differential equation 

(6.5) p{d^)a + p{-d,)b = . 
The following central result holds: 

Theorem 6.3. The map lo = {1^,1^) defined by 

lo-.V^J^, /o(0,?/,0,?/) :=(/_,/+) , f- = a', f+ = b'. 

is one-to-one and extends to an unitary map 

I -n^ Graph(5p) . 

Proof. It is easy to check that the map '■ ^ T is iniective. It 
is surjective too. Indeed, if (/+,/_) G then /o(0, 0, y) = {f+, f-) 
where 

(0,2/,0,|/) = (0,M-it;(0),0,M-ii;(0)), 
0(x) = a{x) + b{—x) , 0(x) = a'(x) + b'{—x) , 

:= / f-{y)dy, b{x):= / f+{y)dy. 



Let us now show that 

II (0, 0, y) III = ||/-(0, 0, y) ||2 + ||/+(0, y, 0, y) ||2 

Since 

and S'p is unitary we need to show that 



||(0,l/,0,^)||| = 2||/-(0,2/,0,y)||' 



2 • 
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Let {(pit, x) , y(t) , (j){t, x) , y(t)) the solution of 

( 4){t,x) - (j)"{t,x) = 
(6.6) <^ y{t)-Ly{t)-w<j)'{t,0^) = 

( e</)'(t,O+) + 0(t,O+)-(w,y(t)) = O 

with initial data {(j){x),y, (l){x),y). Then, by using the evolution equa- 
tion for y and the boundary conditions, one has 

^g^^ j^{\\ym'-{yit),Ly{t))-e\^'{t,o^)\') 

=0*(t, O+)0'(t, 0+) + O+)0(t, 0+) 

and so 



<p*{t, o+)0'(t, 0+) + {<!)')* {t, o+)0(t, 0+) ) dt 

= -{\\yr-{y,Ly)-e\<p'{0^)\') . 
By conservation of energy one has 
2E{(j),y,<i),y) 

' \l + U'\\l + \\yr-{y,Ly)-e\<j>\o^)\' 

'<p*{x)<p{x) + {(f)'y{x)(j)'{x)) dx 



0*(t,O+)0'(t,O+) + (0')*(t,O+)0(t,O+)j dt. 

Inserting in the last equation 4>{x,t) = a(t + x) + h{t — x) and using 
y = Spa', we have 

Ei(l),y,<P,y) = 2\\b'\\l = 2\\I;{(l),y,^,y: 



D 

We define now the maps Rf : V — > 5o(M) by 

:= 4^(0, 2/, 0, 2/) (-a;) 

and the orthogonal spaces Ti^ as the closure, with respect to the energy 
norm, of 

V^:={{(f),y,<i),y)eV : Rf{(f),y,<i),y)eS^{m)}, 

where 

S^iR) := {/ e SoiR) : fix) = 0, ±x < 0} . 

The next theorem show that the subspaces ?i~ and Ti^ are incoming 
and outgoing in the sense of Lax-Phillips scattering theory (see jllj . 
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[T5] . section XI. 11). The proof is a straightforward consequence of the 
previous definitions. 

Theorem 6.4. The subspace H" is incoming and the subspace TC~^ is 
outgoing, i.e. 

u'H- c u'n- , s<t<o, Pi u'n- = {o} , |J um~ = n , 
u'n^ c u'H^ , t > s > , Pi u'n^ = {o} , |J um+ = n . 

t>o teM. 

By the previous theorem and by [Ll\, chapter II, sections 2 and 3, 
there follows 

Theorem 6.5. The unitary maps 

- n^ l2(r) , 

defined as the closures of the maps Rf , provide incoming and outgoing 
translation representations of W , i.e. 

R^U\R^)~^ = T\ 
R^U^ = L2(M±), 

s; = R+{R-)-\ 

where 

T' : L\R) ^ L\R) , T'f{x):=f{x~t). 

Proof. The thesis follows from the Theorem 16.31 and from simple com- 
putations. Otherwise one can use, as we said, the general theory con- 
tained in [Hj. □ 

7. The Lax-Phillips semigroup 

We are now interested in the evolution of the states which are neither 
incoming in the past nor outgoing in the future. To this end one defines 

:= Pf/*P, 

where P is the orthogonal projection onto 

}C:=nein- ®n+). 

Since Ti"*" and ?i~ are orthogonal it is known (see [11], [15] section 
XI. 11) that is a strongly continuous semigroup of contractions on 
/C for positive times: 

Vt>0, Z*:/C^/C, ||Z*||<1, limZ* = 0. 



LAX-PHILLIPS SCATTERING FOR GENERALIZED LAMB MODELS 



27 



The next theorem completely characterizes such a semigroup (let us 
remark that, by fl6.ip . all the roots of p have positive real part): 

Theorem 7.1. The vector space /C is finite dimensional, 

dim/C = deg(p) . 
It is generated by the vectors 



(7.1) 



(i? 



= 0,1, 



- 1. J 



, m . 



rj-r^Q^j^ , for X < 

for a; > ' 

where Zi, . . . are the roots of the polynomial p and z/i,...!/^ the 
respective multiplicities. Moreover 

Z* = e-*^, a{B) = {z^,...,zJi , 

and the matrix representing B with respect to the basis ( [7. ip is the 



Bj where 


5, 


is the 


Uj X 




A. 1 










\ 


z, 


2 




















\0 







- 

h 


1 

/ 



Proof. Since R+ : H ^ L\R) is unitary, R^H^ = L^{R±) and R+{R-)-^ 
S* we have 

By our hypotheses p{iC) 7^ for any ( in the upper complex plane C+. 
Thus by Paley- Wiener theorems the analytic extension to C+ of the 
Fourier transform of f E Sq H S*Sq 

p(-0 



/(C) 



piK) 



has no poles and has zeroes of order z/j at izj, i.e. 



Thus 



C=izj 



2tx y_ 



x^e^'''f{x) = 0, 



= 0,1, 



, i/j- - 1 



L2(M_)nS';L2(M_) = A; = 0, 1, . . . , i/^ - 1, j = l,...,m}" 
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and the finite dimensional subspace K, is whose generated by vectors 
f l7.ip . These vectors are indipendent and so the dimension of K, is 
Yl^=i — degp. To determine the action of on /C is enough to 
note that the evolution of vectors (I7.ip in the outgoing representation 
is given by 



(X - t) = (X - t) V^(^-*)X(-oo,0](x - t) 



Thus 



d 



>kj 



[x-t] 



- {k(t)k-ij{x) + Zj(t)kj{x)) 



and the proof is done. 



□ 
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